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Abstract 

Extrapolating from a sample of 39 AGNs, we examine the impact on the total number of high 
energy Vn induced events (PeV-EeV) expected in IceCube (a lKm 3 "neutrino telescope"), due 
to variations in different parameters involved in the neutrino flux, such as the emission region 
geometry, the estimation models and distributions of the Doppler factor and the variability time. 
This work has been done taking into account different limits of the extragalactic neutrino flux. 

Among our conclusions, we find, in the case of the largest variability time, that the cylindrical 
geometry hypothesis for the emission region, produce a separation of 3a in the total number of 
events relative to the spherical hypothesis. In addition, for similar choices of the burst time, 
spherical geometry and for the upper neutrino flux bound, we obtain a separation of 2.5a in 
the total number of events, for some of the Doppler factor estimations. These differences are 
undistinguishable for other input values. 
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I. INTRODUCTION 



Nowadays, there is a special interest in studying different cosmic rays sources, for instance, 
Active Galactic Nuclei (AGN) [1-3], Gamma Ray Bursts (GRB) [4-7], Dark Matter and 
others [8]. They all help us to understand the structure and formation of our Universe. 
In particular, AGNs are very attractive to investigate, given that they are one of the most 
powerful continuous sources of the highest energy particles in the universe. 

Among these particles, AGNs produce copiously neutrinos. Due to their weak interaction 
with matter, neutrinos keep their flux unaltered in their travel from the AGN to Earth. This 
feature makes them an ideal tool for exploring physical processes within the AGNs. 

At present, our knowledge about the parameters involved in the calculation of the extra- 
galactic neutrino flux is not so accurate. In fact, there is not an absolute agreement about 
some aspects related to these parameters. This is the case of the model used to estimate the 
Doppler boosting factor (e.g. Equipartition [9, 10], Inverse Compton [9, 10] and Variability 
[11] models) and its distribution (power-law-like), the shape of the emission region [12] (e.g. 
spherical and cylindrical) and the variability time scale, which can vary from few minutes 
to many days. A good determination of the Doppler factor is a key ingredient to establish 
the intrinsic photon luminosity, directly related with the neutrino flux. In the same way, the 
shape of the emission region changes the volume considered in the particle density of target 
photons that interact with protons. In addition, the time of burst is inversely proportional 
to the neutrino flux, so that, if the variability time decreases, the neutrino flux increases. 

It is straightforward to note that changes in these parameters affect directly the neutrino 
flux and, hence, imply modifications in the expected value of the total number of neutrinos 
that will be detected on Earth. Consequently, it would be interesting to know the range of 
variation of these events when different assumptions are made about the parameters. If the 
number of events changes in a noticeable way, these numbers can be used as a tool to test 
the model parameters. 

In this paper, we study the impact of these variations in the total number of muon 
neutrino induced events that could be detected in IceCube [13]. This detector is one of 
the so called "neutrino telescopes" (e.g. Amanda, Antares, Nestor) and consists of arrays 
of photomultipliers, located in a Cherenkov medium (ice), under construction in the South 
Pole with a standard volume of lKm 3 . We choose ten years of exposure in the energy range 



2 



of PeV-EeV. In spite of the fact, that this interval must contain a small number of events, it 
is convenient to be studied, considering the absence of background in this energy region. We 
model the neutrino flux based on a catalogue of 39 AGNs and extrapolate for a maximum 
flux given for different proposed limits, which are the Waxman-Bahcall [14], the benchmark 
used by the IceCube Collaboration [15] and the 7-ray bound [16]. 

II. AGNS AND THE UNIFIED MODEL 

Active Galactic Nuclei are one of the brightest known objects in the Universe, with an 
observed luminosity ranging between 10 42 — 10 49 erg/s. According to the Unified Model of 
AGN [17], they have a central engine composed by an accretion disk orbiting a rotating 
super-massive black hole (10 6 to 10 10 solar masses). The AGN is powered by an accretion 
rate of matter onto the super-massive hole of a few solar masses per year, possibly supplied 
by the rotational energy of the black hole. 

The Unified Model also claims other characteristics of the AGNs, such as: two jets 1 
moving in opposite directions and outwards to the accretion disk, an optically thick torus 
and radio lobes. Furthermore, for this model Seyferts, Quasars, Blazars and Radio Galaxies, 
are all AGN, only with different orientations as viewed from Earth, accretion rates and black 
hole masses. In this work, we assume this model including in our test sample all these types 
of galaxies. 

III. NEUTRINO FLUX FROM AGN 

We will consider as our working hypothesis [18] that within inner regions of AGN jets, 
protons are highly accelerated by sheet-like shock waves (blobs, knots or bursts) to energies 
up to the range of PeV-EeV, according to Fermi's stochastic acceleration. These protons 
should then interact with ambient photons emitted by the accretion disk from the thermal 
blackbody radiation (UV bump or blue bump, between 7 x 10 14 Hz to 3 x 10 15 Hz), leading to 
the A + (1232) resonance and followed by the decays (pion from photo-production) nir + — > 
— > V^Vf,. Thus, through this cascade, a neutrino flux is produced. 

1 Generic AGN models do not necessarily involve jets. 
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In order to calculate the observed neutrino flux produced in each AGN, we begin with 
the following generic useful relation 



-EV, 



I E v ^dE v = , (1) 

where the integration limits depend on the range of energies chosen to observe (in our case 
PeV-EeV), L v is the observed neutrino luminosity and di is the luminosity distance to the 
AGN, defined as 

d L = d m (l + z) , (2) 

with z the redshift of the AGN and d m the proper motion distance (see Eq. (B2) in App B). 
The energy spectrum, will be assumed to be proportional to 

§roc^ex P (-^). (3) 

The typical power law energy dependence of the extragalactic neutrino flux is a resemblance 
of the proton spectrum, which is supposed to be accelerated by the Fermi's shock mechanism. 
The slope of the power law, p, is usually taken to be 2, the generic E~ 2 . However, it can be 
considered in the range 1.3—2.7. The spectrum also includes an exponential cutoff depending 
on the AGNs maximum attainable energy, E^ N (we will assume a general cutoff of 10 10 
GeV). 

Therefore, the exact neutrino spectrum formula is 

d ^-MC>^M-m-) dK T E ^(-^) ■ (4) 

We focus our attention on the observable quantities that compose the neutrino luminosity. 
Then, to unfold this luminosity in terms of these observables, we start from the definition 
of luminosity 

L. = ^> , (5, 

^obs 

where N u is the number of produced neutrinos, (E v ) the mean neutrino energy and At f, s 
the observed variability time of the emission region. Now, we give N u from the neutrino 
production sequence explained at the beginning of the section 

N v = (&r A ^„ 7r +)(6r 7r +^jA A , (6) 

where "6r" stands for the branching ratio of the desired decay channel. This gives approxi- 
mately N v a: N w + ~ |Aa- 



It is useful to express the neutrino energy in terms of the energy of the parent proton by 

(E v ) = (x p ^+) (x w ^) (E p ) , (7) 

where (x) represents the average fraction of energy transferred from the parent particle to 
the product, which numerically gives E v ~ \E jX fa \{^E P ). 

On the other hand, to find the number of produced A's we use 

N A = %N y a p ^ A At obs , (8) 

with $ p the initial flux of highly accelerated protons, iV 7 the number of ambient photons 
from the UV bump and a p7T ^ A =10~ 32 m 2 the delta photo-production cross section. $ p and 
N-y are given by 

L P R L^At^ 
P (E P ) V ' ^ (E 7 ) , {9) 
where L p and L 7 are the proton and photon observed luminosities, respectively, (E y ) is the 
mean photon energy (in the present work we take it to be 10 eV), R is the observed size 2 of 
the emission region and V the observed volume of the emission region. It must be noted that 
the size of the emission region and, consequently, its volume depend on four parameters: the 
Doppler factor (5), the variability time scale (At o6s ), the redshift (z) and the angle between 
the velocity vector of the knot and the line of sight (8). The detailed description of R and 
the volume can be found in Sec. Ill A 2. 

Replacing Eqs. (6 - 9) into Eq. (5), we obtain 

_ (br 7T +^ l/ti )( y br A ^ n7T +) (x p ^ 7T +) (^x n ^ v ^ L p L 1 At b s Ro' p ^A , , 

Lv ~ {E^V • (10) 

To write Eq. (10) in a more compact notation, we use the optical depth, r, defined as 

R JjLyAtobs 

T - = K n 7 <T P7 ^A = ti <J P1 -,A , 1-L1J 

where A P7 _,a is the proton interaction length, R, the observed size, which should be under- 
stood as the linear path traveled by the proton in the direction of movement through the 
ambient photons and we write K as 

K = (6r ff +_ H/M )(6r A _„ w +) (x p ^+) (x n ^) « 0.024 , (12) 



In the case of a spherical or cylindrical emission region it corresponds to its radius. 
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where we have used &ta-»7t = 0.5, br\^ u = 0.9998, (x p ^ n ) = 0.2 and (x n ^ v ) = |. 

Since r is a Lorentz invariant, we can go from the observer's frame to the frame attached 
to the knot, so that we can use a simple geometry (see Sec. Ill A 2) to describe the emission 
region, leaving r as follows 

,_ R'(5,At obs ,z,9) L'^At' R'(5,At obs ,z,6) L 7 At o6s 

T T V>(5 : At obs: z,9) (E>) ap ^ A V'(5,At obs ,z,6) (E,) ap ^ A - [i6) 

where primed quantities refer to a quantity measured in the jet and recalling that L ^ At °^ bs is 
itself a Lorentz invariant. This now gives for Eq. (10) 

L u = K r L p . (14) 

However, not all L p will be observed because part of it would have already interacted with 
the ambient photons to produce the A's. Therefore, the intrinsic luminosity will be e T times 
the observed luminosity, L p = e T L Pobs . In our analysis we take L Pobg , which is still unknown, 
as 10% of the photon total luminosity [18]. In addition, we must take into account the 
possible absorption of tc + into the source when the medium is optically thick and therefore 
will not give time for the n + to decay into neutrinos. In this case, we use the substitution 
suggested by Ref. [18], changing r by (1 — e" T ). Thus, the factor (1 — e~ T )e^~ e T ) accounts 
for the absorption of p and n + in the source. Finally the v p and differential flux emitted 
by an individual AGN source can be obtained, using Eq. (4), with 

d ®u u +u u 1 I^Zp ( E v 



dE » /;r;^^exp(-^^) " p v e%»j (15) 

= F(5,V,At obs ,z,L w 6,Er x ,Er n ) E^exp^-J^j , 

where the factor 2 takes into account the contributions of and z/ M 3 . The function F 
includes the parameters that we are going to study 5, V and At obs (the calculation of 6 is 
indirectly affected by the estimation model of S). 



3 The muon antineutrino produced after the fi + decay has almost the same average fraction of energy 
transferred from the n + , {x n ^, Vii ) s=a ix-n^v^), because each lepton shares approximately equal energy. 
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A. Flux parameters to be analyzed 



In this subsection we describe the different hypothesis in the neutrino flux parameters 
studied in this paper. 

1. Doppler boosting factor 

As we have shown in the section before, the Doppler factor, 5, is a piece in the flux formula. 
At present, there are various models for computing 5, which, in some cases, give different 
output values of 5 for the same object. This disagreement could also imply differences in 
the flux that would be interesting to examine. For our work, we consider the following 
estimation models of 5: 

• Equipartition 5 (EQ) [9, 10]: This is calculated assuming the equipartition of energy 
between radiating particles and magnetic field and is obtained from the ratio of the 
observed and the maximum intrinsic brightness temperatures. The maximum intrinsic 
temperature is supposed to be the equipartition temperature. 

• Inverse Compton 5 (IC) [9, 10]: In this case, the Doppler factor is determined from the 
comparison of the observed and predicted X-ray fluxes, assuming that this is caused 
by inverse Compton scattering of synchrotron photons off the radiating particles. 

• Variability 5 (Var) [11]: This is obtained in a similar way to 5eq, but is weaker than 
this, since it takes the third root of temperatures' ratio. Additionally, it replaces 
the observed brightness temperature instead of the variability brightness temperature, 
which is estimated from a total flux density flares associated with VLBI components 
emerging from the AGN core, covering nearly 20 years of variations. 

• Minimum 5 (Min) [19]: This model is computed only from the proper motion data 
using the minimum 5 value allowed by 



where T is the Lorentz factor and (3 app the apparent transverse velocity in units of 
c (the speed of light). We are going to use this model independent estimation as a 
reference. 




(16) 
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Due to different proper motion observations of underlying shocks for a single jet, there 
should be a rather broad range of Doppler factors [20]. Then, it follows that a certain 
distribution in 5 must be taken into account when computing the neutrino flux. According 
to Ref . [20] it is consistent with a steep power law distribution of intrinsic Doppler factors. In 
our calculation, we fit with a power law (5~ 7 ) the data obtained for each type of calculation 
of S, as described in Appendix. B. In Sec. V we explain how the neutrino flux is averaged 
over the fitted distribution. 

2. Emission Region Geometry 

In the study of the emission region geometry, we consider two possibilities for its shape 
in the proper frame: a spherical and a cylindrical volume. These models are based on the 
observed variability time At obs [12]. However, we should state that At obs could give, at 
most, the estimation of only one of the dimension of this region and therefore the other two 
should be inferred somehow from this. In this sense, the formulas given here are actually 
not relativistic transformations (as in App. A), but rather estimations and should not be 
misunderstood. 

To understand the next formulas, it is appropriate to define a quantity D = that 
includes both the Doppler factor (5) and the cosmological effects given by the redshift (z). 
The geometries are described as follows: 

• Spherical geometry: it has two options for the jet-frame source radius R! . The first 
option, that we call SphereA, uses the common formula 



The second option, that we call by SphereB, uses a refined formula for R', that 
consider the angle to the line of sight, 9, and the shock velocity, (3 shock (used value 
/3 s /iocfc=0.5). This expression is given by [12] 



R' = DcAt obs . 



(17) 



R' = 



DcAt, 



obs 



(18) 




TD(co8(6) - /3) + 



Dsm(9) 

r 
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Cylindrical geometry: with radius R! and jet-frame length I. The region is rapidly 
energized by a plane shock traveling along the jet, such that photons are emitted 
immediately after shock passage from a thin disk-like region downstream of the shock. 
From here on, we assume I = aR', with a variable (used value a = 2). 

R' = 



DcAt obs 



f Abs 



P.. 



TD(cos(9) - /3) +Dsin(9) 



(19) 



3. Variability Time Scale 



The observed fluctuations over the whole spectrum in the AGNs, or variability time scales 
of the bursts, are one of the main characteristics of Active Galaxies [21]. This variability 
suggested that the emitting regions were extremely compact, leading to the hypothesis that 
AGN were powered by massive black holes. 

The time scale varies according to the class of AGN and waveband of the spectrum. The 
shortest fluctuation observed in an AGN is attributed to Mrk421 with a 15 minutes duration. 
However, there have been variability on time scales of months to years [21]. In the case of 
blazars, a rapid variability or flaring is common. Observed short time scales of a day or more 
in the X-ray band are characteristic to TeV blazars (e.g. Mrk421, Mrk501, PKS2155-304). 
They have even substructures with shorter time scales of 10 3 — 10 5 sec [22]. On the other 
hand, many EGRET sources had time scales of several days [18] and the optical waveband 
flares are fairly slow. 

After this considerations, we choose the variability time scales ranging from 15 minutes 
to 10 days. 



IV. OUR AGN SAMPLE 



Although there are more than 876 Blazars, 11777 Seyfert-1 and 48921 Quasi Stellar 
Objects catalogued by Veron-Cetty (2003) [23], only few of them contain the parameters 
important for our interests. After a thorough search of AGN catalogs, we have built our 
own AGN sample from some of them, so that gathering their data, we included relevant 
information to our investigation, such as the proper motion of jet components, luminosity 
and different Doppler factor estimations. From the combination of Ref.[20], the latest and 
widest AGN survey available on proper motion observations, with Ref.[24], which includes 
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luminosities, and Ref.[9, 11], containing Doppler factor estimations, we have selected only 
the objects that have information on all parameters, remaining a sample of 39 AGN. In 
order to get a proper motion distribution for the AGN's jets we add more data in the proper 
motion observations from earlier surveys [10, 19, 25-28]. 

We must note that these multiple observations are not from the same epochs, which might 
introduce some uncertainties. However, we think that it is good enough for our purposes. 

A. Description of the sample 

Our total sample of 39 AGN, is divided into tree sub-classes: 10 BL Lacs (blazars) (see 
Table I), 27 Core-dominated Quasars (see Table II) and 2 Radio Galaxies (see Table III). 
Tables [I, II, III] show the observational data that describes the properties of each source (i.e. 
AGN model parameters). The data is presented as follows: 

• Column(l) International Astronomical Unit source designation [Name B1950]. 

• Column(2) Alternative source Name. 

• Column(3) Redshift. 

• Column(4) Monochromatic luminosity at 15 GHz, from [24], completed with [29]. 

• Column(5) Component identifier of the knot. 

• Column(6) Angular radial speed: proper motion. 

• Column(7) Equipartition Doppler factor, from [9]. 

• Column(8) Inverse Compton Doppler factor, from [9]. 

• Column(9) Variability Doppler factor, from [11]. 

• Column(lO) Reference for proper motion. 

A comment is in order, in our calculations of the neutrino flux, we need to use the UV 
luminosity L lobs in Eq. (15) (see Sec. III). For that reason we extrapolate the luminosity 
given in this sample, which is in the radio component of the spectrum, to find the UV 
luminosity. This estimate has been done in the following way: from the luminosity fractions 
(for different parts of the electromagnetic spectrum) of the set of AGNs presented in Ref . [30] , 
we compute the average luminosity fractions, such that the radio part represents 0.52% and 
the optical-UV region 26%. Then, with the relation L uv as 50L rad we do the extrapolation. 
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B. Sample characteristics 



To better understand the characteristics of the AGN sample, we present in this subsection 
several histograms and comparisons of its most relevant properties. 

We display in Fig. 1 the angular distribution, in right ascension and declination, of our 
sample's 39 AGNs. Most of the sources are located in the northern hemisphere due to 
the abundance of observatories in that region. Nevertheless, this bias does not affect our 
purposes, since we are going to consider the neutrino flux as isotropic. 

We show in Fig. 2 the extrapolated UV luminosity (in the observer's frame) distribution 
of our sample. It can be seen that Quasars have the highest luminosity, with a mean value of 
7.23 x 10 46 erg/sec, while Blazars and Radio- Galaxies have 4.82 x 10 45 erg/sec and 4.50 x 10 45 
erg/sec, respectively. These observations are consistent with our expectations, since Quasars 
are thought to be the most luminous objects in the universe. 

In Fig. 3 we present the extrapolated UV luminosity versus its intrinsic counterpart, 
which was obtained using Eq. (A3), for each Doppler factor estimation model. From that 
equation we expect that in the case of a 5 slightly greater than one, the observed luminosity 
should be higher than the intrinsic luminosity. This situation is reflected in our sample, 
where the observed luminosities for each AGN is higher than their corresponding intrinsic 
values, which is represented by the points below the straight line. This indicates that most 
AGN jets in the sample are moving towards Earth and suffer the Doppler boosting. Only 
few AGN, in the case of the Equipartition and Inverse Compton estimation models, have 
the opposite behavior. This is because, for these models, some Doppler factor estimates are 
lower than one. On the contrary, for the Variability model all Doppler factor estimates are 
greater than one, meaning that their luminosities will suffer the so called Doppler boosting. 

The apparent transverse velocity distribution for each knot in the sample is shown in 
Fig. 4. From this plot it is clear to note that Quasars have higher values of apparent 
transverse velocities than Blazars and Radio Galaxies, which is in agreement with Fig. 2, 
due to the relationship between Doppler factor, which depends on the apparent transverse 
velocity, and luminosity. We also notice here that some of the AGN jets are moving away 
from Earth, which match with our observations in Fig. 3. 

The distribution of the mean angle to the line of sight, calculated using the method 
described in Appendix B, is given in Fig. 5 for different Doppler factor estimation models. 



11 



As we could have supposed from the previous plots, we see that the majority of the AGN jets 
are pointing almost directly towards us. This fact does not entirely reflects the assumptions 
of the Unified Model, which predicts, in average, greater angles for Quasars and Radio- 
Galaxies, than we have in our histograms. Meanwhile, Blazars are compatible with small 
angles, which is expected, since their jets are pointing almost directly to us. The mean 
9 values found for the whole sample are for 9 E q ~ 5.1°, 9 IC ~ 9.4°, 9 Var « 7.6°, and 

9rnin ~ 13.4°. 

In Fig. 6 we plot the distribution of the Doppler factor for each estimation model. We note 
that the Equipartition and Inverse Compton cases give lower 5 values than the Variability 
case. Their mean values are < 5eq >= 8.00, < 5jc >= 7.65 and < 5v ar >— 9.61 and they 
range between 5 EQ = [0.14,29], 5 IC = [0.15,27] and 5 Var = [0.98,26.21]. We found that 
these histograms can be fitted with the expected power law distribution. After applying the 
distribution discussed in Appendix B we obtain: < Seq >= 8.11, < 5[c >= 6.54, < by ar >= 
8.49 and < 5 m i n >= 7.04. These new values correspond to a power law distribution with 
negative slope (5~ 7 ) with exponent ^eq = 0.38, 77c = 0.51, ^var = 0.53, and 7 m j n = 0.55. 

Finally, we compare in Fig. 7, using two-dimensional plots, the different types of esti- 
mations for the Doppler factor, 5eq, he an d 5y ar . We find that Seq and Sjc are strongly 
correlated, in contrast, <5y ar is uncorrected with any of them. 



V. ^-INDUCED EVENT NUMBER CALCULATION 

We calculate the number of AGN-neutrino induced events observed on Earth within the 
context of the proposed IceCube [13] detector. As we are interested in astrophysical point 
sources (i.e. AGNs), we use muon neutrinos because they can achieve the highest angular 
resolution (0.7°) due to the large muon track, which allows to reconstruct their direction 
[31]. Moreover, //^-induced events have a higher effective volume [31] than other neutrino 
flavors. We restrict our search to the highest energy events in the PeV-TeV range, within 
IceCube's sensitivity [13], thus we avoid the background signals from atmospheric neutrinos. 
Then, the number of + induced muon events is calculated by [8] 

N, = T dQ dE Vii ( (E^) ) A eff P surv (E Vll , 9) P v ^ 9, E$*) , 

J JlQiS eV \ I tot 

(20) 
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where T is the exposure time (we use 10 years in our calculations) and A e ff is the effective 
area of the detector, taken from Ref. [15]. Since, for this energy range, IceCube is capable 
of observing both hemispheres, we integrate over the whole celestial sphere. Notice that in 
this expression, corresponds to the zenith angle (9 = is pointing to the South Pole) and 
should not be confused with the angle to the line of sight used to describe the neutrino flux 
in the previous sections. 

The total isotropic flux of AGN neutrinos summed over all sources is defined as 

= <?»£ / d6 ^3,mod F 3 (6, V, At, Zj, L, ob3j , (9) {j)mod , 10 18 eV, 10 15 ev) 

j=l \ °minj 

= Cn(0)^x^" Pex p(-^v) > 

(21) 

where {<f>) tat represents the sum of the terms Fj (see Eq. (15)), after averaged over the 
probability density function (pdf) in 5, defined for each j source of the total 39 AGN sources 
of our sample. The probability density function varies for each estimation model (mod) and 
is given by 

pdf(5) j>mod = - 7 , (22) 

f, 3 dS 

where 7 corresponds to each Doppler factor estimation model (^eq — 0.49, = 0.58, 
Ivar = 0.55, 7 m j„ = 0.76). The integration limits, S maX j and S min j, are chosen for a single 
AGN source, taking the maximum (S maX j) and minimum (5 m inj) of its corresponding set 
4 of Doppler factors (these limits depend on the 5 estimation model). The normalization 
constant C n in Eq. (21) allows us to extrapolate the total flux of our small sample to the 
maximum possible. This constant is given by 

f~i Abound (O'i') 

1 T r max i \^"- > ) 
Wtot 

where (pbound stands for the three different limits that we have chosen: the Waxman-Bahcall 
limit 3 x 10~ 8 GeV^cm^s^sr -1 [14], the benchmark used by the IceCube Collab. 1 x 10~ 7 
GeV~ 1 cm _2 s~ 1 sr- 1 [15] and the 7-ray bound 6 x 10~ 7 GeV~ 1 cm _2 s~ 1 sr- 1 [16]. We must 



4 When there is only one value for an AGN source, we assume an average range. 
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note that the given bounds are still below the limits set by recent experiments and above 
IceCube's sensitivity [15]. The limit for energies between 1 PeV and 3 EeV from the latest 
AMANDA results [32] derived with a 90% C.L. is E 2 $ V {E) < 0.99 x 1(T 6 GeV^cm" V^r" 1 
and the AMANDA-II-2000 data is still being processed. ((fyToT a ^ 1 GeV, is the maximum 
value among the set of values of {4>) tot obtained for all the variations performed in our work 
(different types of S estimation, class of geometry emission region and variability time of the 
burst). It is relevant to mention that there is a possibility that these sources were optically 
thick enough to absorb most of the protons, only allowing neutrinos to scape, which would 
indeed exceed the conservative WB bound. 

Continuing with the description of Eq. (20) we define P sur v (E Vll ,9), the survival proba- 
bility of a neutrino crossing the Earth, by 



X(0) 



(24) 



Psurv ( E »^°) = ex P 

where x (@) is the column density given by 

x (0)= / p(r(9,l))dl, (25) 
Jo 

with, L , the neutrino distance traveled through Earth calculated using 



Z„„ (0) = \/ {Re ~ d) 2 cos 2 (0) + 2dR E - d 2 - (R E - d) cos (9) , (26) 

where R E is the Earth's radius (6378 Km) and d the distance to the detector (for IceCube 
it is 1.4 Km). 

In Eq. (25) we integrate over the neutrino path crossing the Earth using the parameterized 
Earth's density, p (r), given by the Preliminary Earth Model (PREM [33]) and the inclusive 
interaction length, L*°* t , is given by 

^ " o% ( E \) N Av ' (27) 

where the inclusive cross section, a*^, is taken from Ref. [34], which includes the charged 
and neutral current contributions. 

The probability that a muon produced in a charged current interaction reaches the 
detector before losing all its energy down to a value lower than the detector's energy threshold 
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E thr (10 3 GeV for IceCube) is defined by 5 



<?cc (E Vft ) 



1 



1 — exp 




where the differential charged current cross section, j s taken from Ref. [35] and the muon 
range, R^, is taken from Ref. [36]. 

We include the effects of neutrino oscillations [37] by multiplying Eq. (20) by 0.5, given 
that the ratio of the initial to the oscillated muon neutrino flux, 0^ : 0°* c , would be 2:1. 

It is worth mentioning that there are other scenarios where muons can lose a significant 
amount of energy before their decay, equivalent to considering only a single [38] , otherwise 
we would have to consider half of the neutrino flux (without neutrino oscillations), which is 
different from the one we assume here. 

VI. RESULTS 

We present in Fig. 8, Fig. 9 and Fig. 10, an analysis of the dependence of the expected 
total number of induced and events versus the variability time (ranging from 15 minutes 
to 10 days) for different Doppler factor estimation models, emission region geometries, and 
for ten years of exposure in IceCube. We use the generic spectrum E~ 2 and normalize the 
figures with the WB limit (Fig. 8), the IceCube Collab. benchmark (Fig. 9) and the 7-ray 
bound (Fig. 10). In each figure we show four plots related to the different type of 5, where 
we display four curves corresponding to the following hypotheses in the geometry of the 
emission region: SphereB, Cylinder, SphereB-D and Cylinder-D (the last two are averaged 
over the 5 distribution). 

In general, we observe that the tendencies of the three figures are the same and they only 
differ in the scale of number of events. The behaviour in At b s of the curves within each 
figure, can be understood from the factor (1 — e~ T )e^~ e T ) in Eq. (15), with r oc l/At obs . 
Also, using this factor, we can comprehend two important features of our plots. First, we 
see why all the graphs converge to the same number of events (given for the flux limit) for 
an small At b s , which implies that for this time we are not able to distinguish among the 

5 For down-going muons if R^{0) > d, then we chose = d, because the Earth's density traversed by 

neutrinos is smaller. 
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different model parameters under study. Second, we observe the hierarchy of the curves, 
where the number of events is higher for the cylinder than for the sphere geometry. This is 
not strange since the spherical hypothesis gives a greater volume than the cylindrical one, 
which implies fewer number of events given that r oc 1/V. Furthermore, comparing the 
number of events for the same geometry, we obtain slightly higher values in the case where 
we average over the Doppler factor distribution; this is because the distribution includes 
lower values of S, which gives, in general, a greater number of events. 

In order to quantify the differences observed in our plots, at the largest variability time, in 
the number of events, for the different parameter hypotheses studied (e.g. type of geometries 
and Doppler factor calculations), we have defined the following estimator: 



which gives us these discrepancies in terms of the number of cr's, with Nf?™* the total 
number of expected neutrinos, where limit stands for the bound used in the extragalactic 
neutrino flux and hyp(i) refers to the various hypotheses compared. 

Our first analysis has the purpose to evaluate how large the differences in the number of 
events are for the spherical and cylindrical geometry. This is shown in Table IV, where we 
have substituted in Eq. (29) hyp^ = SphereB(SphereB-D) and hyp^ = Cylinder (Cylinder- 
D), evaluated for each type of Doppler calculation. From this Table, we can see that, in 
the case of the WB limit, the differences go from 0.75a to 1.53(7, when we do not include 
the 5 distribution. Meanwhile, when the S distribution is considered, these differences are 
even smaller going from 0.53(7 to 0.89er. For the IceCube benchmark, the differences are 
significant and range between 1.36-2. 79a and 1.02-1.62(7 for the cases without and with S 
distribution, respectively. For the 7-ray bound the differences are greater, ranging from 
3.34-6.84cr and 2.35-3.98(7 as described before. 

In addition, we have done an analysis similar to the one before, only that in this case 
we have calculated the differences in the number of events considering the six different pairs 
of Doppler calculation's types. This is displayed in Table V and Table VI, for a specific 
geometry, spherical and cylindrical, respectively. There are only two favourable situations 
for discriminating the Doppler factor models using the number of events. This is for the case 
of a spherical geometry, the 7-ray bound and without 5 distribution, where we have obtained 



hyp(i) h VP(j) 



(29) 
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2.67a in separation for the events predicted using Seq and Sic- The other case is when we 
compare 5ic and 5v ar , where we have got 3.05er. For all the rest of pair comparisons we 
have got negligible differences, less than 1.6a in the splitting of the number of events. 

In spite of not including variations in "a" , the relation between the radius and the length 
in the cylinder (I = ra), we have checked that these changes will not lead to a relevant 
difference in the number of events. 

As complementary information, we show in Fig. 11 the total number of expected events, 
for the three flux bounds, as a function of the slope p of the AGN spectrum, ranging from 
1.3 to 2.7. It is clear that in p = 2 the number of events in this figure coincides with the 
maximum number of events in our precedent figures. Thus, from this plot we can rescale 
the results given in our tables, only multiplying the N a by N ^(p) / N ^(p = 2), which at 
p = 2.7 is around \/2. 

VII. CONCLUSIONS 

We have studied different neutrino production scenarios in AGNs for the highest energy 
neutrinos (PeV-Eev) and their impact in the total expected z/^-induced events in the IceCube 
detector, due to variations in the emission region geometry, the Doppler factor estimation 
model and the variability time. We have obtained the total neutrino flux based on a sample 
of 39 AGN (i.e. point sources observations), where the limits in the total neutrino flux were 
the WB bound, the IceCube Collab. benchmark and the 7-ray bound. 

We have found as a global behaviour that if most AGN had long times of burst, which 
seems natural to assume, we should be able to distinguish, in some cases, among the different 
models. Otherwise, for small times of burst, no information about these model parameters 
could be extracted through neutrino observation in a detector with the capabilities of Ice- 
Cube. 

More specifically, for a 10-day fluctuation time and the 7-ray flux limit (6 x 1CT 7 
GeV _1 cm _2 s _1 sr _1 ), we could determine the emission region geometry after 10 years of 
IceCube exposure. We observe a separation greater than 3<r in the total number of events, 
between the cylindrical and spherical geometry, for all the 5 cases. On the other hand, if the 
flux bound was below the WB limit (3 x 10~ 8 GeV _1 cm _2 s~ 1 sr _1 ), the geometry would not 
introduce major changes in the total number of events. When the benchmark is 1 x 10~ 7 
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GeV _1 cm -2 s _1 sr _1 , we are able to discriminate between geometries in half of the 5 cases, 
without including the 5 distribution, with a separation in the total number of events greater 
than 2(7. 

It is important to note that it is possible to infere from Fig. 8, Fig. 9 and Fig. 10, any 
intermediate case (a mixture of both geometries), between the cylindrical and spherical 
geometry for the emission region. In this sense, how far or close we are from some of these 
extreme behaviours (purely cylindrical or spherical) is a measure of the tendency of the 
sample. 

For the Doppler factor estimation models, in most cases, the total number of neutrino 
induced events gives no relevant information or, in other words, the differences in 5 do not 
affect this number. Nonetheless, for a spherical emission region, At f, s =10 days, the 7- 
ray bound and without considering the 5 distribution, it is possible to distinguish between 
8eq-8ic an d Sjc-Svar obtaining an event separation greater than 2.5a. 

It is relevant to keep on mind that our analysis has considered that the neutrino spectrum 
had slope p = 2. Therefore, for a slope p > 2, the separation in terms of N a must increase, 
which means that we are going to enlarge the number of cases where we can differentiate 
the geometries and the Doppler factor estimations models. 

As closing remarks, we must say that our work establish possible tendencies in the pa- 
rameters in the AGNs. In fact, we expect that these tendencies become more robust and 
measurable in the future, when we have a bigger sample containing the necessary information 
and accurate point sources observations. 
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APPENDIX A: RELATIVISTIC TRANSFORMATIONS 



For the energy and time transformation between rest-frame to observer-frame, which in- 
cludes not only the common relativistic transformation but also the cosmological expansion, 
we use 

E' = { ^-E oha , (Al) 

At' = j—jAt obs , (A2) 

where 5 is the relativistic Doppler factor, also known as bulk Lorentz factor and z the AGN 
redshift 6 . 

There are luminosity transformations [17] which consider the observed integrated spec- 
trum and solid angle distributions for moving, isotropic sources and for a continuous standard 
conical jets. However, in order to proceed in a more general way (i.e. model independent) 
we consider the following jet photon luminosity transformation 

L' j= {^jhobs. (A3) 

APPENDIX B: OUTFLOW ANGLE, DOPPLER FACTORS AND DISTRIBU- 
TIONS 

In order to obtain a more accurate number in the exponent, 7, in the power law formula 
(<5~ 7 ) for each Doppler Factor estimation model distribution, we enlarge our 39 data of the 
Doppler factor estimation model (there is one for each source) calculating the Doppler factor 
for each knot in the source. The method for calculating these Doppler factors is described 
in the following lines. 

First, we calculate the apparent transverse velocity P app (j, k) corresponding to each knot 
(index k) for its AGN source (index j), given by 7 

^ U , k) = !^M, (Bl) 



6 There should not be any misleading impression with the fact that the redshift means a receding object 
and most Doppler factors values an approaching source, because the redshift is due to cosmological effects 
and the Doppler factor to the moving knots in the jets. 

7 The formula from Pearson&Zensus (1987) widely used in older papers is replaced by Eqs. (B1,B2) which 
consider a non-zero cosmological constant. 
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where k) is the proper motion for each knot in the source, d m (j) is the "proper motion 
distance" to the source, calculated using [39] 

d m (j) = 77" / V(l + zY{l + n m z) - z(2 + z)n A dz , (B2) 
-"0 Jo 

with the Hubble constant H = 71 Km S' 1 Mpc' 1 , the matter density VL m = 0.27 ± 0.04, 
and dark energy density 1)a = 0.73 ± 0.04, according to WMAP latest results [40]. 

Then, we calculate the angle to the line of sight, 0(j,k) mod , for each proper motion 
component and its corresponding Doppler factor for each estimation model (the index mod 
= EQ, IC, Var), using 



k) mo d — arctan . 
In the case of Minimum Doppler factor, which gives the minimum angle, we use 



(B3) 



0(j, k) min = arccot [P apP (j, k)] . (B4) 

^From the distribution of 9(j,k) mo d for each source, we obtain its mean value {0)^ mod 
and assign it to the source, which is well justified [24]. Thus, we can construct the Doppler 
factor distribution for each model in terms of the Doppler factor of the knots for each AGN 
source. 

6V, k) rnod = - n f~^ k) Z d , (B5) 

1- P(j,k) mod cos (9) {j)mod 

where 

P{J ' k)mod = WlWW^ • (B6) 

We discard Doppler factors that could correspond to jets affected by significant bends, 
twists and inward moving features. 
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APPENDIX C: THE AGN SAMPLE: MODEL PARAMETERS 

TABLE I: BL Lacs 
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TABLE I: - Continued 
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TABLE II: Quasars 
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TABLE II: - Continued 



Source Alias Redshift Lum. Comp. /i r Seq Sic <Var Ref. 

(W Hz' 1 ) (mas yr' 1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

B 0.15 V94 

0336-019 CTA 26 0.852 3 x 10 27 C4 0.18 13 12 19.01 J01 

C3 0.42 J01 

B 0.22 K04 

C 0.21 K04 

0420-014 OA 129 0.915 6.3 x 10 27 B 0.20 21 13 11.72 J01 

B 0.03 K04 

C 0.29 K04 

0528+134 PKS 0528+134 2.07 4.6 x 10 28 K2/U2 0.16 0.71 2 14.22 H00 

F2 0.40 J01 

B1+F3 0.26 J01 

B2+F4 0.33 J01 

B3 0.19 J01 

B4 0.15 J01 

B 0.08 K04 

0552+398 DA 193 2.365 3.6 x 10 28 . . . 0.06 1.2 2.2 14.2 G93 

0.06 G97 

0.04 V94 

0804+499 1.432 1.9 x 10 27 B 0.13 21 16 26.21 K04 

0836+710 4C 71.07 2.17 1.39 x 10 28 . . . 0.35 8.5 6.7 10.67 G93 

0.27 G97 
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B 0.23 V94 

D 0.14 V94 
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. . . 1.15 V94 

1226+023 3C 273 0.158 1.5 x 10 27 C3 0.79 8.3 4.6 5.71 C88 

C4 0.99 C88 

C5 1.20 C88 

C7a 0.76 C88 

1.62 G93 

0.88 G97 

K10/U10 0.77 H00 

K9/U9 0.94 H00 

K8/U8 1.15 H00 

K7/U7 1.06 H00 

K4/U4 0.99 H00 

B5+E2 0.33 J01 

B4+E1 0.85 J01 

B3 0.66 J01 

B2 0.61 J01 

Bl+Gl 0.68 J01 

H5+K6+D1 1.60 J01 
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TABLE II: - Continued 
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(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 








H3+K5+D2 


0.70 








J01 








B 


1.05 








K04 








c 


1.36 








K04 








D 


0.83 








K04 








E 


1.27 








K04 








F 


0.79 








K04 








G 


0.41 








K04 








I 


0.88 








K04 








C2 


1.15 








V94 








C3 


0.79 








V94 








C4 


0.99 








V94 








C5 


1.20 








V94 








C7 


0.65 








V94 








C7a 


0.76 








V94 








C8 


0.92 








V94 








C9 


0.82 








V94 


1253-055 


3C 279 


0.538 1.3 x 10 28 




0.50 


14 


14 


16.77 


C88 








B2 


0.11 








C88 










0.68 








G93 










0.14 








G97 








K4/U4 


0.17 








H00 








Kl/Ul 


0.25 








H00 








B3 


0.18 








J01 








B2 


0.17 








J01 








E2+B1 


0.27 








J01 








D 


0.31 








J01 








B 


0.28 








K04 








C3 


0.12 








V94 










0.50 








V94 


1308+326 


OR 017 


0.997 5.8 x 10 27 


K2/U2 


0.19 


4.3 


5.2 


11.38 


H00 








Kl/Ul 


0.96 








H00 










0.45 








G97 








Kl 


0.13 








V94 








K2 


0.75 








V94 








K3 


0.29 








V94 








B 


0.31 








K04 


1510-089 


4C 09.56 


0.36 3.4 x 10 26 


K3/U3 


0.45 


10 


11 


13.18 


H00 








Bl 


0.51 








J01 








Dl 


0.28 








J01 








D2 


0.63 








J01 








B 


0.85 








K04 








c 


0.57 








K04 


1633+382 


4C 38.41 


1.807 6.5 x 10 27 


B3 


0.14 


0.83 


2.2 


8.83 


J01 








Bl 


0.20 








J01 








B 


0.15 








K04 








Q 


0.10 








K04 


1641+399 


3C 345 


0.594 5.9 x 10 27 


C2 


0.48 


1.5 


4.1 


7.45 


C88 








C3 


0.30 








C88 








C4 


0.30 








C88 










0.66 








G93 










0.33 








G97 








B 


0.49 








K04 








C 


0.37 








K04 








C2 


0.47 








V94 








C3 


0.30 








V94 
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TABLE II: - Continued 



Source 


Alias 


Redshift Lum. 


Comp. 


Mr 


I7V 1 




S v 


Rcf. 






(W Hz- 1 ) 




(mas vr~^) 










(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 








C4 


0.31 








V94 








C5 


0.23 








V94 


1928+738 


4C 73.18 


0.303 6.2 x 10 26 


Al-4 


0.60 


3.4 


3.4 


3.71 


C88 










0.81 








G93 










0.39 








G97 








B 


0.29 








K04 








c 


0.30 








K04 








E 


0.12 








K04 








Al 


0.32 








V94 








B 


0.37 








V94 








c 


0.34 








V94 








C2 


0.51 








V94 








C3 


0.57 








V94 








C4 


0.40 








V94 








C6 


0.40 








V94 








C7 


0.60 








V94 








C9 


0.31 








V94 


2021+614 


OW 637 


0.227 2.6 x 10 26 




0.04 


0.9 


1.1 


1.59 


C88 










0.05 








G93 










0.02 








G97 








B 


0.04 








K04 








c 


0.05 








K04 










0.02 








V94 


2134+004 


PHL 61 


1.932 2.9 x 10 28 




0.01 


15 


27 


11.49 


C88 










0.01 








G93 








B 


0.02 








K04 


2145+067 


4C 06.69 


0.999 1.1 x 10 28 


B 


-0.01 


15 


21 


7.81 


K04 








c 


0.03 








K04 








D 


0.03 








K04 


2223-052 


3C 446 


1.404 1.2 x 10 28 




0.14 


18 


16 


11.39 


G93 








B 


0.49 








K04 








c 


0.31 








K04 


2230+114 


CTA 102 


1.037 4.3 x 10 27 




0.65 


1 


1.5 


14.23 


C88 










0.69 








G93 








B3 


0.25 








J01 








B2 


0.34 








J01 








Bl 


0.33 








J01 








B 


0.03 








K04 








c 


-0.05 








K04 








E 


-0.23 








K04 


2251+158 


3C 454.3 


0.859 1.2 x 10 28 


2 


0.05 


5.6 


4.6 


21.84 


C88 








4 


0.35 








C88 










0.48 








G93 










0.09 








G97 








B3 


0.14 








J01 








B2 


0.34 








J01 








Bl 


0.53 








J01 








B 


0.04 








K04 








2 


0.05 








V94 








3 


-0.05 








V94 








4 


0.35 








V94 








5 


0.21 








V94 
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TABLE III: Radio Galaxy & Seyfert 



Source 


Alias 


Rcdshift 


Lum. 


Comp. 


fir 


<5bQ 


Sic 


Svar 


Ref. 








(W Hz^ 1 ) 




(mas yr ) 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1845+797 


3C 390.3 


0.057 


2.9 x 10 24 




0.98 


0.37 


0.38 


1.16 


G93 










B 


0.54 








K04 










C 


0.60 








K04 


0430+052 


3C 120 


0.033 


8.2 x 10 24 


A 


1.35 


ll 


4.1 


0.98 


C88 










B 


2.53 








C88 










C 


2.47 








C88 










D 


2.66 








C88 










E 


2.54 








C88 












3.55 








G93 












2.06 








G97 










K1B/U1B 


1.62 








H00 










T S" -I A /TT1 A 

KIA/UIA 


2.22 








H00 










B 


1.77 








K04 










C 


1.80 








K04 










D 


1.36 








K04 










G 


1.59 








K04 










I 


l.5l 








K04 










H 


2.08 








K04 










A 


1.35 








V94 










B 


2.53 








V94 










C 


2.47 








V94 










D 


2.66 








V94 










E 


2.54 








V94 



Note.- The acronyms used in the references for the proper motion, given in column (10), correspond to: 
C88 ([27]), G93 ([19]), G97 ([10]), H00 ([25]), J01 ([26]), K04 ([20]), and V94 ([28]). 
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FIG. 1: Distribution of the AGNs sample given in right ascension and declination angles. 
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FIG. 2: Distribution of the observed UV-Luminosity in the sample. 
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FIG. 3: Intrinsic versus observed luminosities according to the Doppler estimation model. 
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FIG. 4: Apparent transverse velocity distribution in the sample. 
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FIG. 5: Distribution of the mean angle to the line of sight in the sample. 
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FIG. 6: Doppler factor model distribution in the sample. 




FIG. 7: Comparison of the different types of Doppler factors for the sample. 
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FIG. 8: Expected total AGN z^-induced events as a function of the variability time, for differ- 
ent assumptions in the model parameters studied. The maximum flux is normalized with the 
Waxman-Bahcall limit [14] and p = 2 is taken in the energy spectrum. "D" stands for the use of 
5 distribution. 
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FIG. 9: Expected total AGN z^-induced events as a function of the variability time, for different 
assumptions in the model parameters studied. The maximum flux is normalized with the bench- 
mark used by the IceCube Collab. [15] and p = 2 is taken in the energy spectrum. "£>" stands for 
the use of 5 distribution. 
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FIG. 10: Expected total AGN ^-induced events as a function of the variability time, for different 
assumptions in the model parameters studied. The maximum flux is normalized with the 7-ray 
bound [16] and p = 2 is taken in the energy spectrum. "D" stands for the use of 5 distribution. 
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v flux limit 


Without 8 distribution 

8eq $ic 8y ar Smih 


Using 8 distribution 

$EQ Sic Svar SMin 


WB 


1.53 0.75 1.45 1.08 


0.81 0.56 0.89 0.53 


IceCube 


2.79 1.36 2.64 1.97 


1.48 1.02 1.62 0.96 


7 — ray 


6.84 3.34 6.47 4.82 


3.62 2.51 3.98 2.35 



TABLE IV: Separation in terms of N a between the number of events obtained for the spherical and 
cylindrical geometries, using the three v flux limits, all 8 cases and At b s = lOdays. 



v flux limit 






Without 8 


distribution 






EQ-IC 


EQ-Var 


EQ-Min 


IC-Var IC-Min 


Var-Min 


WB 


0.68 


0.04 


0.35 


0.59 


0.31 


0.31 


IceCube 


1.25 


0.08 


0.64 


1.09 


0.56 


0.56 


7 — ray 


3.05 


0.18 


1.57 


2.67 


1.38 


1.38 




Using 8 distribution 


WB 


0.33 


0.08 


0.26 


0.24 


0.07 


0.18 


IceCube 


0.61 


0.15 


0.48 


0.44 


0.13 


0.32 


7 — ray 


1.49 


0.36 


1.17 


1.09 


0.32 


0.79 



TABLE V: Separation in terms of N a between the number of events obtained for different pairs of 8 cases 
and the three v flux limits, using the spherical geometry and At b s = lOdays. 



v flux limit 






Without 8 


distribution 






EQ-IC 


EQ-Var 


EQ-Min 


IC-Var IC-Min 


Var-Min 


WB 


0.04 


0.03 


0.05 


0.01 


0.01 


0.02 


IceCube 


0.07 


0.05 


0.09 


0.02 


0.02 


0.04 


7 — ray 


0.18 


0.13 


0.23 


0.04 


0.05 


0.09 




Using 8 distribution 


WB 


0.09 


0.16 


0.01 


0.06 


0.10 


0.16 


IceCube 


0.17 


0.29 


0.02 


0.11 


0.19 


0.30 


7 — ray 


0.43 


0.70 


0.04 


0.27 


0.47 


0.73 



TABLE VI: Separation in terms of N a between the number of events obtained for different pairs of 8 cases 
and the three v flux limits, using the cylindrical geometry and At b s = lOdays. 
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FIG. 11: Expected total AGN muon neutrino induced events normalized with different bounds of 
the v flux limit as a function of the exponent in the power law assumed in the AGN-z^ flux. 
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